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Abstract 

We study pair annihilations of the neutrahno dark matter in the minimal supersymmetric 
standard model with CP violation. We consider the case that the higgsino mass and the 
trilinear scalar couplings have CP-violating phases of order unity, taking a scenario that 
the scalar fermions in the first two generations are much heavier than those in the third 
generation to avoid a severe constraint from experimental limits on electric dipole moments. 
It is found that, when the lightest neutralino (x) is bino-like, the cross sections of the 
VF-boson pair production XX ~^ W^W~ and the lightest Higgs boson pair production XX 
— > H2H2 for nonrelativistic neutralinos can be significantly enhanced by the phase of the 
higgsino mass. The relic density of the neutralino can be considerably suppressed by this 
effect. However, even this suppression is not enough to make bino-like dark matter consistent 
with a cosmological constraint. We also discuss the effect of CP violation on the positron 
flux from neutralino pair annihilations in the galactic halo. 
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I. INTRODUCTION 



Existence of a considerable amount of cold dark matter (CDM) in the present 
Universe has been regarded as a robust ingredient in recent astrophysics and cos- 
mology Q|. In particular, the Wilkinson Microwave Anisotropy Probe (WMAP) has 
determined the relic abundance of CDM as 

Ocdm/^' = 0.1126+Ho9, (1) 

where ^cdm is the CDM energy density normalized by the critical density and h ^ 0.7 
is the Hubble constant in units of 100 km/sec/Mpc 3]. The relic density of the CDM 
is expected to be determined with more accuracy by analyses of increasing WMAP 
data and forthcoming data from the future project Planck j^. 

Identification of the CDM as a particle still remains to be settled. Among possible 
candidates, the hghtest superparticle (LSP) in supersymmetric models is one of the 
most attractive candidates In supersymmetric models, conservation of a discrete 

symmetry called R-parity is nessesary to avoid rapid proton decay, but it in turn plays 
a crucial role to guarantee the stability of the LSP. In the minimal supersymmetric 
standard model (MSSM) f?!, the LSP is typically the lightest neutralino \o\ which is 
a linear combination of neutral gauginos and higgsinos 

X = Xi = NnB + Ni2W' + N^^H^, + N^Hl (2) 

where B is the U{1)y gaugino (bino), is the neutral SU{2)l gaugino (wino), and 
and H2 are the two neutral higgsinos with opposite hypercharges. The coefficients 
Nii {i = 1,2,3,4) are the elements of the 4x4 unitary matrix which diagonalizes 
the neutralino mass matrix 

A number of theoretical analyses on the relic density of the neutralino CDM in the 
MSSM have already been done extensively BHll,|l2,13,14,15,16,17,18,Q- 
However, most analyses have been done with an assumption of no CP violation in 
supersymmetric parameters. CP violation is not only observed in particle physics but 
also expected to play an essential role to explain baryon asymmetry in the Universe. 



Since it is known that the standard source of CP violation in the Cabibbo-Kobayashi- 
Maskawa mixing matrix is not enough to produce required asymmetry of order 10"^'', 
supersymmetric CP violation is expected to gerenate the correct amount of baryon 
asymmetry 21 1. Thus, it is generically important to include the effect of super- 
symmetric CP violation. When the supersymmetric CP phases are of order unity, they 
typically lead to too large electric dipole moments of the neutron, the electron and 
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Hg atom to satisfy experimental limits 



24| if masses of superparticles are 



within a TeV scale. However, if superparticles in the first two generations are much 
heavier than 1 TeV, the supersymmetric CP phases of order unity are still allowed. 

On the other hand, a complete analysis of the neutralino relic density f^^^^ includ- 
ing all the relevant effects in the MSSM with CP violation is still missing. The effect 
of CP violation on the neutralino pair annihilation cross section into the fermion pairs 
through the sfermion exchange was studied in Ref. 

In Ref. Q, 

the effects of a 

CP-violaing phase in the trilinear scalar couplings on the neutralino pair annihilation 
cross section were examined, including all the final states and taking into account 
the scalar and pseudoscalar mixing in the neutral Higgs sector js?, 28|. Recently, the 
dependence of Q^h"^ on CP-violating phases through supersymmetric loop corrections 
to the bottom-quark mass was examined in Ref. 29 1 . Partial wave treatment in the 
presence of CP violation was studied in Ref. |2Ql]. An analysis related with electroweak 
baryogenesis was done in Ref. l8ll . 

n 

In this paper, we extend our previous analysis and perform a reanalysis on 
the effects of CP-violating complex phases in the MSSM on pair annihilation cross 
sections of the lightest neutralino and its relic density in the present Universe. We 
consider the case that the higgsino mass and the trilinear scalar couplings have CP- 
violating phases of order unity, taking a scenario that the scalar fermions in the first 
two generations are much heavier than those in the third generation to avoid a severe 
constraint from experimental limits on electric dipole moments. We include all the 
contributions to the neutralino annihilation cross section at the tree level, taking into 
account the scalar and pseudoscalar mixing in the neutral Higgs sector. In the absence 
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of CP violation, it is known that fermion pair productions XX ~^ ff usually give the 
largest contribution to the total cross section for a bino-like LSP. In the present paper, 
we shall show that, unlike the case without CP violation, the W^-boson pair production 
XX ~^ W^W~ or the lightest Higgs boson pair production xX -^^2-^2 ^^'^ S^^^ the 
largest contribution in the presence of CP violation. We also discuss the effect of CP 
violation on the positron flux from neutralino pair annihilations in the galactic halo. 

This paper is organized as follows. In section m we describe the structure of the 
MSSM with supersymmetric CP violation, and present the relevant interactions. In 
section IIIH we discuss the effect of CP violation on the cross section of the lightest 
neutralino pair annihilation. In section llVj we briefly review the formalism to compute 
the relic density of the neutralino. In section|3 we summarize the procedure to obtain 
the positron flux. In section IVIl we present our numerical results. Finally concluding 
remarks are given in Section FVIII 

II. THE MSSM WITH SUPERSYMMETRIC CP VIOLATION 

In the present work, we consider a general MSSM [3] in which all the interactions 
are specified by the following input parameters at the weak scale 

Ml, M2, M3, /i, mo. A, tan/3, m^, (3) 

where Mi, M2 and M3 are the mass parameters for the bino, the wino and the gluino, 
respectively. The parameter fi represents the Higgs mixing mass, and tuq is the super- 
symmetry breaking common mass parameter for the sfermions in the third generation. 
The corresponding mass parameters for the first two generations are assumed to be 
10 TeV to suppress the electric dipole moments ^. The parameter A is the common 
trilinear scalar coupling for the third generation, while the ones for the first two gen- 
erations are neglected. The ratio of the vacuum expectation values of the two neutral 

^ The sfermion mass spectrum considered here naturahy arises in a minimal supersymmetric 5*0(10) 
GUT model 0. 
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Higgs fields is denoted by tan/?. The quantity rriA is a parameter which coincides 

|. In eq. (0), 



with the pseudoscalar Higgs mass in the absence of CP violation [27 
Mi {i = 1,2,3), fi and A can have CP-violating phases in general. For the gaugino 
masses, we assume the GUT relation 

3Mi ^ M2 ^ Ms , . 

5 g'' g^ gr ^ ^ 

where g', g and gs are the gauge coupling constants for f/(l)y, SU{2)l and SU{?>)c 
gauge groups, respectively. With this relation, the neutralino LSP can be bino-like 
(I Ml I -C higgsino-like (|Mi| > or their mixture (|Mi| ^ 

In the present analysis, we assume that the gaugino masses Mj are real, and ex- 
amine the effects of the CP-violating phases of /i and A 



H = \fi\exp{i9^), A = \A\exp{i9A] 



(5) 



These phases induce imaginary parts in the mass matrices for the neutralinos, the 
charginos and the sfermions. They also induce mixings between scalar Higgs fields 
and a pseudoscalar Higgs field through radiative corrections 0, Q. The MSSM 
contains two scalar neutral Higgses 0? ^i^d 0° with hypercharges = —^(02) 

= — I and a pseudoscalar neutral Higgs ^4° as physical fields. In general, nonzero 
phases in /i or A induce (p^-A^ and 0°^^° mixings at the one-loop level, so that the 
3x3 mass-squared matrix A4'jj for the neutral Higgs fields has to be diagonalized. 
We calculate the Higgs mass-squared matrix with the one-loop effective potential, 
including the contributions of the third generation fermions and sfermions. The mass 
eigenstates (r = 1, 2, 3) are related with the CP eigenstates (05,02,^°) by a 3x3 



rotation matrix O as follows 



-"2 



/ 0?^ 



o 



(6) 



The Higgs boson mass eigenvalues are obtained as OA4^fjO^ = diag(m^o,'^^o,''^^o)- 
In eq. (jH)), H2 is defined as the lightest Higgs. The heavier Higgses and are 



defined such that |0i3| < IO33I. In the case of 9^ = 9a = 0, it follows that {H2) is 
the heavier (lighter) scalar Higgs field, and iJg is the pseudoscalar Higgs field. Note 
that the parameter itla in eq. (jHl) is not a mass eigenvalue in general. The scalar- 
pseudoscalar mixing in the neutral Higgs sector implies O^s 7^ {r = 1,2) so that 
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interactions of the neutral Higgs bosons are significantly modified 

Our phase convention in the Higgs sector is defined as follows. In addition to 9^ and 
9ai the coefficient of a Higgs mixing term (0?02) can have a complex phase 9i2 (See 
Ref. j27[). In general these phases induce a nonvanishing relative phase ^ between 
the vacuum expectation values of the two Higgs fields through the relevant tadpole 
minimum condition. However, 9i2 and ^ are not separately physical quantities, and 
only their sum ^ + 9i2 is rephasing invariant. We adopt a convention that ^ = at 
the one-loop level so that nonvanishing 6*12 will be induced by the tadpole condition 
in the presence of nonvanishing 9^ or 9a- 

In the present work, we focus on the VT-boson pair production XX ~^ W~^W~ 
and the lightest Higgs boson pair production xX ~^ H2H2 . For later discussion, we 
explicitly describe the structure of the interactions relevant for these processes in the 
following. 

When the Higgsino mass fi has a nonvanishing phase, the elements of the neutralino 
mixing matrix have an imaginary part in general. The interactions of the neutral 
Higgses with the neutralinos x^ = 2, 3, 4) in the presence of CP violation have 
the following structure 

= 5 E E X? {C^"^"' - Cf^"'-,^ X?/^?. (7) 

r=l i,j=l 

In the absence of CP violation, either the scalar coupling Cs or the pseudoscalar 
coupling Cp in eq. ((Zj) is vanishing for every mass eigenstate Namely, the 

scalar coupling C5 is vanishing for H^, and the pseudoscalar coupling Cp is vanishing 
for and //^ • ^^e presence of CP violation, however, both Cs and Cp are nonzero 
for every mass eigenstate. These couplings can be written as 

Cf"' " = Re(G^^), 



Cp'''' ^ = zIm(G;^), (8) 



where 



= \{g'N,^-gN,2)[N,:,Ori-N^,Or2 

+ 20,3 (Ar,3 sin/5 - cos/5)] + (i ^ j). (9) 

Note that the couphngs in eq. (jHl) are vanishing if the hghtest neutrahno x (= Xi) is 
a pure bino or a pure higgsino. However, such a case can not take place in practice 
since there always exists a finite mixing between the bino and the higgsino. 
The VT-boson-W^-boson-neutral Higgs boson interaction is given by 

£^^^0 = ^ (10) 

r=l 

The expression for the coefficient is 

(jWWH^^ = ^mvF[0,icos/5 + 0,2sin/3], (11) 

where is a mass of the W^-boson. Numerically, the magnitude of (7^^^° is the 
largest for the lightest Higgs if^ > the others are much smaller as we explicitly see 
in Section IVTl 

The chargino-neutralino-iy-boson interactions are given by 

i^..-w. = {cI^'^'^'Pl + C^^'^'^'Pn) XW; + h.c. , (12) 

k=l 

where Pl = {I - 75)/2, Pr = {1 + 75)/2 and 



of = g [sJIn,,v:, - N,,V:^ . (13) 



The 2x2 unitary matrices U and V diagonalize the chargino mass matrix M.y- as 



U*M.^~V ^ = diag(m^-, m^-) [l^. Note that these couphngs are vanishing in a pure 
bino limit. 
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Finally, trilinear couplings for the neutral Higgs fields are given by 



(14) 



r,s,t=l 



The coupling constant C*^°^°^° is written as 



4 cos 9w 



(Orl COS P - Or2 slu /?) (0,20*2 - O.iOn + OM3 COS 2/?) 



+ (the other permutations of r, s, t) 



(15) 



where mz is the Z-boson mass, and 6w denotes the Weinberg angle. 

III. NEUTRALINO PAIR ANNIHILATIONS WITH CP VIOLATION 

Neutralino pair annihilations are of importance when we examine the relic density 
of the neutralino and its indirect detection. The complete expressions for the total 
cross section of neutralino pair annihilations in the case of no CP violation are found 

Ref. Q. We have extended analvs. of Ref. Q to inco.po.ate *e effee. of 
CP violation, and derived the full expressions for the neutralino annihilation cross 
section at the tree level with CP-violating phases, taking into account the modified 
interactions including (jZj), (fTn|l and (fT^ . Calculation of the total pair annihilation 
cross section a involves a number of final states: 



XX 



ff, iJX(^,s = l,2,3), H^H- 



W^W, ZZ, W^H^, ZH^.{r = 1,2,3). 



(16) 



The complete set of the analytic expressions including all the final states are quite 
lengthy and complicated, so we do not list the full result. Instead, we present only 
some of the results for the W^-boson pair production xx ~^ W'^W~ and the lightest 
Higgs boson pair production xX ~^ -^^2-^2 which are essential in later discussion. In 
our numerical analysis in section IVH however, we include all the contributions to 
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the neutralino annihilation cross section at the tree leveL Note that we can assume 
nonrelativistic neutrahnos in the analysis of positron flux, while we have to keep the 
exact form to calculate the relic density. 

Before presenting the expressions with CP violation, let us briefly summarize some 
generic features on the cross sections in the case without CP violation. When the LSP 
is bino-like, the fermion pair production XX ~^ ff (/ = u, c, t, ■ ■ ■) typically gives a 
dominant contribution to the total pair annihilation cross section for nonrelativistic 
neutralinos. There are three diagrams which contribute to this process: neutral Higgs 
boson exchange, Z-boson exchange and sfermion exchange (see Fig. QJ. The s-wave 
amplitude of this process for nonrelativistic neutralinos is always suppressed by a 
fermion mass (except for the top quark). Because of the s-wave suppression, the cross 
section of a light fermion pair production from nonrelativistic neutralino annihilations 
is much smaller than that for a heavy fermion. If the neutralino is heavier than 
the ly-boson and/or the lightest Higgs boson, the processes xX ~^ W~^W~ and/or 
XX ~^ H2H2 open up. However, for a bino-like LSP, these processes give smaller 
contributions than that for the // final state in the case without CP violation. 

On the other hand, when the LSP is higgsino-like, the W-hoson pair production 
XX ~^ W^W~ gives a dominant contribution if the neutralino is heavy enough. This 
process involves three diagrams: neutral Higgs boson exchange, Z-boson exchange and 
chargino exchange (see Fig.|21). Since this contribution has no s-wave suppression, the 
neutralino pair annihilation cross section is much enhanced so that the resultant relic 
density of the neutralino is typically too small to satisfy the WMAP constraint (^. 

However, in the presence of CP violation, new features appear in the W^-boson pair 
production for a bino-like LSP. In order to see this fact, let us present essential parts 
of the analytic results for the cross section of xX ~^ W^W~ in the following. 

The analytic expression for the neutral Higgs exchange contribution to xX WW 
is given by 

'^WW'^ ~ on 2 o 4 
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X 



(s — 4m, 



3 qWWHOqXXH° 
r=l 



Pho(s) 



+ s 



r=l 



c 



WWH°(jXXH° 



PHsis) 



;i7) 



where 



Pho (s) = s- m]jo + i Tho niHO , 



(18) 



niy is a mass of the hghtest neutrahno, Tjjo denotes the decay width of H^, P' 



w 



a/I — 4m^/s is the velocity of the ly-boson in the center of mass frame, and s is 
a Mandelstam variable. The relative velocity between the two colliding neutralinos 
is denoted by v. Numerically, the contribution from the lightest Higgs exchange (r 
= 2) is dominant as expected from eq. (fTT|) . Note that C^^'^ becomes nonzero in the 
presence of CP violation. For nonrelativistic neutralinos s — > 4m^, the first term in the 
square bracket of eq. (jl7|) vanishes. On the other hand, the second term includes the 
factor s rather than s—Am?. Therefore this term does not vanish in the nonrelativistic 

A. 

limit once CP violation is turned on. 

Similarly, the chargino exchange contribution and the interference term between the 
Higgs and the chargino exchange diagrams for XX ~^ W^W~ contain contributions 
which do not vanish for s 4m^ in the presence of CP violation. There are also 
other contributions to the 14^-boson pair production which include Z-boson exchange 
and possible interference terms, but we neglect them for the moment since they are 
typically subdominant. With all the relevant contibutions together, the expression 
for the total cross section of xX ~^ W^W~ for nonrelativistic neutralinos {v 0) is 
given by 

2 



awwv 



Svr 



{ml - mlr] 



k=l 



A 



w 



+ 



4m^ — im'^ym'^ + 3m^ 



w 



2m^ 



E 



r=l 



(JWWH'^qXX.H? 

Pm{^ml) 



" Xk 



E 

k=l 



(19) 
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where 



m 



w 



mr 



m? _ and 

Xfc 



C 



i Im 



L 



xw- 



xtxw- 

R 



(jxtxw- 



(20) 



In eq. fll9|) . the second term in the square bracket represents the effect of CP viola- 
tion, while the first term does not vanish even without CP violation. Note that we 
include Z-boson exchange and possible interference terms in the numerical analysis 
even though we neglected them in eq. ()19j] . 

Likewise, the lightest Higgs boson pair production xx ~^ -^^2-^2 involves contribu- 
tions which are drastically enhanced by CP violation. This process occurs via two 
diagrams: neutral Higgs boson exchange and neutralino exchange (see Fig. El). The 
expression for the total cross section of this process for nonrelativistic neutralinos {v 
^ 0) is given by 



where A,. ^ 



mr 



Phi 
IGvr 



E 

r=l 



(JH°H°H0(JXXH? 



2 = 1 



A^ 



(21) 



m^o, Ph^ denotes the velocity of H2 in the center of mass 



frame, m^o is the mass of the neutralino , and 



D ? 



2i Im 



X°XH^2 f(jxUH°2 



(22) 



Note that the cross section (PT|) is vanishing in the absence of CP violation. 

The above expressions (|19p and (|2H) are valid only for the nonrelativistic limit s 
Am?. In our numerical calculation, however, we use the exact expression for general 



s instead of the above expressions. 



IV. RELIC DENSITY OF THE NEUTRALINO 

In this section, we briefly review the formalism to compute the neutralino relic 
density in the MSSM with supersymmetric CP violation jsj. The time evolution 
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of the neutralino number density in the expanding Universe is described by the 
Boltzmann equation 



where H is the Hubble expansion rate, n'^ is the number density which the neutrahno 
would have in thermal equilibrium, and a is the total cross section of the neutralino 
pair annihilation into ordinary particles. The quantity (erf) represents a thermal 
average of av. 

In the early Universe, the neutralino is assumed to be in thermal equilibrium where 
= n'^. As the Universe expands, the neutralino annihilation process freezes out, 

and after that the number of the neutralinos in a comoving volume remains constant. 

Using an approximate solution to eq. (j^ . the relic energy density = m^n^ at 

present is given by 



where x = T/m^ is a temperature of the neutralino normalized by its mass, (?* (~ 81) 
represents the effective number of degrees of freedom at freeze-out, and Gn denotes 
the Newton's constant. and are the present temperatures of the neutralino and 
the photon, respectively. The suppression factor (T^/T^)^ ^ 1/20 follows from the 
entropy conservation in a comoving volume js^. The value of x at freeze-out, xp, is 
obtained by solving the following equation iteratively: 



Typically one finds xp ~ 1/20 which implies that the neutralinos are nonrelativistic 
at freeze-out. 

The thermal average in eq. (PHj) should be carefully treated for accurate calculation 
of the relic density. In literatures, expansion of the thermal average in powers of the 
temperature (av) ^ a + bx is widely used. However, it is known that the expansion 




(23) 




(24) 




(25) 
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breaks down when a varies rapidly with the energy of the neutra. 
gerenal, one has to use the exact expression for the thermal average 

1 



inos, hence, in 



3 



{av) 



MTKlimjT) 



^ ^ dsa{s){s - AmD^sK^ , (26) 



X 

where Ki {i = 1, 2) are the modified Bessel functions. Note that we should not take 
the nonrelativistic limit f to calculate the cross section a in the integrand. In 
our analysis, we perform a numerical evaluation of the exact thermal average ()26|) to 
obtain the relic density accurately. 

In the present analyses, we neglect coannihilation effects 
although they are crucial when the next-lightest superparticle is nearly degenerated 
with the LSP. The investigation of the effect of CP violation on the coannihilation 
cross sections is left for future work. 

V. POSITRONS FROM NEUTRALINO ANNIHILATIONS IN THE 
GALACTIC HALO 

Cosmic ray observations provide interesting probes for indirect detection of CDM. 
Even though the neutralino pair annihilations described in section IIVI have already 
frozen out at present, neutralinos gravitationally accumulated in the galactic halo 
still can pair-annihilate to produce, e.g., cosmic 7-rays, neutrinos, antiprotons and 
positrions 0, 0] . Among various such observations, the High- Energy Antimatter Tele- 
scope (HEAT) reported an excess of cosmic ray positrons Previous analyses on 
the positron flux in the MSSM without CP violation have suggested that some en- 
hancement mechanism such as nontrivial distribution of CDM and/or so-called boost 
actors is necessary in order to explain the excess by neutralino LSP annihilations 

3 44, Q- 



We describe our procedure to evaluate the positrion flux from neutralino pair an- 
nihilations in the galactic halo in the following. The positron flux we measure can be 
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written as 



44| 



dE 



j deG{E,e)J2i^,v)Me), (27) 



where v = 10~^ is a typical neutralino velocity in the galactic halo, and po = 0.43 
GeV/cni'^ represents a local halo dark matter density. Positrons produced in the halo 
are decelerated or accelerated during propagation in the halo until they are detected. 
The energy of a positron at detection is denoted by E, while that at production is 
denoted by e. The function G{E, e) is a Green's function which describes propagation 
of positrons in the galactic halo. In our calculation, we use the Green's function in 
Ref. 2 with a containment time r = 10'^ yr. The cross section cij represents that for 
the process XX ~^ h where the symbol i runs over every possible final state: i = ff, 
W'^W^, etc. Computation of the positron flux requires the cross sections cxj only for 
V = 10~^. In our numerical calculation of cij for v = 10~^, we use s ~ 4m^(l + f ^). 

The function /j(e) is the positron energy spectrum at production which originates 
from the process xX ~^ There are many contributions to /j(e), and some of them 
including hadronization are not calculable. In order to perform a detailed analysis, one 
has to utilize a Monte Carlo simulation code such as PYTHIA Q]. In the present 
analysis, however, we content ourselves with a rough estimation to calculate /i(e). 



employing the methods used in Refs. j43 



44| as follows. 



The most energetic positron comes from the direct production xX ~^ e^e~ where 
the positron has the maximal energy e ^ m^. However the cross section is extremely 
small due to the s-wave suppression by the electron mass so that this contribution 
is invisible. On the other hand, heavy fermion productions xX ~^ cc, bb, t^t^ and 
ti, may produce enough positrons as decay products (e.g., b — > cW^* followed by 
W^* e^z/g), where the typical positron energy is e ~ Also, hadronization 

of quarks results in a shower of charged pions, and the contribution of the positrons 
from charged pion decay (vr"*" — >■ m"*" — > e"*") are evaluated using the data from e'^e~ 



collider experiments as in Ref. 44 1. 



For a neutralino heavier than the W^-boson, positrons with energy e ~ can 
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be produced by the VT-boson production xx ~^ W^W' followed by the decay —>■ 
e+z/g. There is also continuum positron radiation from muons, r-leptons and heavy 
quarks produced in the ly-boson decay (i.e., ^ fi'^ ^ e^, ^ r+ ^ e'^, etc.). 
The contribution of positrons from decay of charged pions produced in the ly-boson 
pair production is estimated as in Ref . j43| . Similarly, contributions from the Z-boson 
pair production xX ~^ ZZ followed by the decays Z — > e"'"e~, etc. are also included. 
It is known that, when the LSP is higgsino-like, the gauge boson final states W~^W~ 



and ZZ can give rise to enough positron excess 



. For the process including H2 in 



the final state, positrons from H2 — * bb are taken into account. 

For the cosmic ray electron flux, we use = 0.07 {E/ GeV)~^'^ 



cm ^sr ^GeV ^sec ^ j47[. The background flux of cosmic ray positrons are estimated 

dF+/dE 
dF+/dE + dF^/dE 



With ,,^f/,:'lf,_^,j, = 0.02 + 0.10(E/GeV)-o-5 



VI. NUMERICAL RESULTS 

In this section, we present our numerical results. In all the figures, we fix some of 
the parameters in eq. (jHI) as 

mo =\A\ = 1 TeV, niA = 500 GeV, tan/3 = 5, (28) 

and consider variations of the remaining parameters as follows 

e^ = o-7r, eA = o-n, 

M2 = 100 - 700 GeV, = 100 - 700 GeV. (29) 

In what follows, we shall show that the process xX ~^ W^W^ or xX ~^ -^^2-^2 can 
give a significant contribution to crv in the presence of CP violation. For this purpose, 
let us first present the magnitude of the relevant couplings in this process. 

We begin with a discussion of the t- and u-channel chargino exchange diagram 
(see Fig. ^ which is known to give a dominant contribution to xX ~* W^W~ in the 
absence of CP violation. In Fig. |3J the quantities and D'^^ defined in eq. ()2()j] 
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involved in the lighter chargino exchange contribution to XX ~^ W^W~ are plotted as 
a function of the CP violating phase 6^ for M2 = 300 GeV, = 400 GeV and 9a = 
with the other parameters fixed as in eq. (j2Hl)- The dashed and solid curves correspond 
to the result for and D^^, respectively. It is seen that both and have 
nontrivial dependence on 6^. The magnitudes of the coefficients in eq. ^ for this 
choice of parameters are obtained as \Nii\^ ^ 0.97, |A''i2p < 0.003, |A''i3p < 0.02 and 
|-^i4p ^ 0.006, hence the LSP is bino-like. Of course, the Higgs exchange diagram is 
absent if x is a pure bino. However, x actually has small but nonvanishing higgsino 
components due to bino-higgsino mixings even when Mi is much smaller than 
The small higgsino components turn out to play a crucial role in the calculation of 
neutralino pair annihilation cross sections. 

In the presence of CP violation, the largest contribution to xX ~^ W^W~ can be 
given by the neutral Higgs exchange diagram which includes three coupling constants 
(jwwH'^^ (jxxH? ^xxH° ^ 2, 3). In Fig. U the dependence of the absolute 
value of the coupling constants C^^^'' on 9^ is presented for the same choice of 
parameters as Fig. |3J It is seen that the coupling C^^^° is the largest for the 
lightest Higgs ifg independent of 

On the other hand, the imaginary parts of the pure imaginary coupling constants 
C^^-'- (r= 1,2, 3) for the same choice of parameters are shown in Fig. El Without CP 
violation, only the coupling for is nonvanishing. In the presence of CP violation, 
however, the couplings for and H2 also become nonvanishing and even comparable 
to that for H^. As these figures clearly show, the lightest Higgs exchange diagram 
gives the dominant contribution among the three Higgs exchange contributions. We 
do not show a result for the coupling , since it is not essential for our discussion. 

In Fig. the cross section times relative velocity av for v = 10~^ is shown for the 
same choice of parameters as Fig. |31 The solid, long dash-dot, short dashed, short 
dash-dot and short dash-long dash lines correspond to the contributions from W^W~, 
ZZ, bb, T^T~ and H2H2 final states, respectively. The bold solid line represents the 
sum of all the contributions. Without CP violation {9^ = 0, tt), a fermionic final state. 
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bb, is dominant as usually expected for a bino-like LSP. In this case, the cross section 
of XX ~^ W^W~ is dominated by the chargino exchange diagram in Fig. |21and the 
Higgs and Z-boson exchange contributions are negligible. The resultant contribution 
of the W^W~ final state is, however, smaller than the bb contribution. For 9^ 
7r/2, the lightest Higgs exchange contribution to the W~^W~ production is much 
enhanced due to the effect of CP violation to dominate over the chargino exchange 
contibution. Because of this enhancement, the W^W~ contribution can be larger 
than the bb contribution. Similar enhancement due to CP violation can be seen for 
ZZ and H^H^ final states. It is found that the W+W' or H^H^ final state gives the 
largest contribution for O.lvr < < O.Qvr. 

In Fig. ISl the relic density ^^^"^ shown as a function of for the same choice 
of parameters as Fig. 0] In the shaded region, the relic density is consistent with the 
2a allowed range of the WMAP result 

0.094 < n^h^ < 0.129. (30) 

It is found that the relic density is considerably suppressed for 9^ ^ 11/2. However, the 
relic density is still too large to satisfy the WMAP constraint even for 9^k, ti/2. In the 
presence of CP violation, the cross section is certainly much enhanced, but even this 
enhancement can not make bino-like dark matter consistent with the cosmological 
constraint. In order to make the bino-like LSP cosmologically allowed, one has to 
resort to some other mechanism such as resonant annihilations or coannihilations by 
tuning some parameter. 

The 9a dependence of av for v = 10^^ is shown in Fig. IHl for the same choice of 
parameters as Fig. |3] but 6'^ = 0. It is seen that all the relevant contributions are 
insensitive to 9a- As a result, the relic density is also insensitive to 9a {^^^"^ ~ 3.3). 
We have not found any significant dependence on 9a even if we vary M2 and ^ for 
our choice of the other parameters in eq. ()28p . 

When the higgsino components are increased, the relic density shows different 
behavior. Fig. [101 represents the 9fj^ dependence of av for v = 10~^ in the case of 
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M2 = 300 GeV and = 200 GeV where the LSP has sizable higgsino components. 
The values of the other parameters are the same as Fig. |3] The magnitudes of the 
coefficients in eq. ((21) for this choice of parameters are obtained as jA'^nP ~ 0.65, 
I^i2p ^ 0.03, |A''i3p ^ 0.2 and |A^i4p ~ 0.1, hence the LSP can be regarded as a 
mixture of the bino and the higgsinos. For such a mixed LSP, the magnitude of the 
couplings C^^ and C^^ in eq. (fT^ are increased compared with the case of a 
bino-like LSP. Because of this enhancement, the W^W~ final state gives a dominant 
contribution av ~ 10~^ GeV~^ for any 0^. The ZZ and H2H2 final states are also 
much enhanced and either of them gives the second largest contribution. However, 
the cross sections for the gauge boson final states W^W~ and ZZ are less sensitive 
to for a mixed LSP than the case of a bino-like LSP, even though the H2H2 final 
state still has an enhancement for 6fj_ ~ 7i/2. Therefore, the total cross section shows 
only mild dependence on 6^. One might expect that the W^W~ final state is also 
enhanced as in the case for the H2H2 final state. However, what actually happens 
is that there occurs a considerable cancellation between the Higgs exhange and the 
chargino exchange diagrams for a mixed LSP ^. Because of this cancellation, the 
W^W~ final state shows only mild dependence on 9^. Similar cancellation occurs for 
the ZZ final state as well. 

The relic density for the same choice of parameters as Fig. ^1 is shown in Fig. El 
For a mixed LSP, the 9^ dependence of f^^^^ much weaker than the case of a bino- 
like LSP. However, it is found that the allowed region appear at 9^ ^ 2tt/3. Thus 
CP violation can be essential to find cosmologically allowed regions in the parameter 
space. 

Let us briefiy comment on a higgsino-like LSP. When x is higgsino-like, the cou- 
plings Cf^ and in eq. (fT^ have a maximal magnitude of order g without 
any other small suppression factor, and almost insensitive to 9^. In this case, the W- 

^ In the first version of this paper, the sign of the interference term between the Higgs exhange and 
the chargino exchange diagrams in the process XX ^ W^W~ was wrong. Hence this cancellation 
was not seen in the first version. 
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boson pair production gives the dominant contribution, and this leads to a much larger 
cross section av ~ 10~® GeV~^ than that for a bino-like or a mixed LSP. Therefore 
the relic density for a higgsino-like LSP is too small to satisfy the WMAP constraint, 
and almost insensitive to 6^. 

Let us present a global map in the M2) plane. In Fig.[T21 the final states giving 
the largest contribution to av for v = 10~^ are displayed in the M2) plane for 
mo = |A| = 1 TeV, uia = 500 GeV, tan/3 = 5, 9a = and 9^ = 0. The regions where 
the final states W+W', bb, ti and W^H^ give the largest contribution are shown in 



different gray scales. T 



le white region is excluded by the LEP limit on the chargino 



mass m- > 104 GeV 



49l | and the lightest Higgs mass ruj^o > 113 GeV 
bino-like LSP (M2 <^ the bb final state typically gives the largest contribution, 
while for a higgsino-like LSP (M2 ^ the W^W~ final state is the largest. When 
the LSP is relatively heavy, the tt final state can be the largest. When both M2 and 
\fi\ are large (~ 700 GeV), the IV^H"^ final state is the largest. In the darkest strip 
in Fig. El the relic density is consistent with the WMAP 2a constraint. In the region 
M2 ~ 500 GeV, where 2m^ ^ f^H^y there occurs resonant annihilation to the bb and 
ti final states via heavy Higgs boson exchange, so that the WMAP allowed region 
is extended to a large region for M2 ~ 500 GeV. It is seen that aside from the 
resonant annihilation region, the cosmological constraint (j30|) is satisfied for a mixed 
LSP. For a bino-like LSP, the relic density is too large to satisfy the cosmological 
constraint ()3()j) . while it is too small for a higgsino-like LSP. 

The similar result to Fig. but for 9^ = tt /2 is shown in Fig. ^1 It is found 
that, unlike the case without CP violation, the W^W~ production is significantly 
enhanced and gives the largest contribution even for a bino-like LSP. Note, however, 
that even this enhancement of the W~^W~ final state can not save the bino-like region, 
as the WMAP allowed region does not appear for M2 < |/i| except for the resonant 
annihilation region in Fig. El Also, there appears a region where the H2H2 final 
state gives the largest contribution. 

In Fig. El variation of the relic density with 9n normalized by that for 9n = 0, 
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^x(^m)/^x(^a« ~ 0)' is shown as a function of M2 for mo = |A| = 1 TeV, = 
500 GeV, tan/3 = 5 and 9a = 0. Varing 9^ in the range < < vr, the rehc density 
hes between the two sohd hnes for |/i| = 200 GeV. The region between the two dashed 
hnes represents the corresponding resuh for = 600 GeV. It is confirmed from this 
figure that the variation of the rehc density with 9^ is typically small for a higgsino-like 
LSP. 

Finally we illustrate the effect of CP violation on the positron flux from neutralino 
annihilations in the galactic halo. The positron fraction e+/(e+ + e~) versus positron 
energy Ee+ is shown in Fig. ^1 for the same choice of parameters as in Fig. E3 
The LSP in this case is a mixture of the bino and the higgsinos. The dashed and 
solid lines correspond to the results for 6*^ = and 9^ = tt/2, respectively. The 
points with error bars represent the data from HEAT measurement ,42]. For both 
9f^ = and 9fj^ = 7r/2, the positron flux is dominated by the background so that 
the supersymmetric contribution from the LSP annihilations is almost invisible. The 
similar result for M2 = 600 GeV with the other parameters fixed as in Fig. ^1 is 
shown in Fig. ^| This choice of parameters corresponds to a higgsino-like LSP. In 
this case, the supersymmetric contribution gives a large excess due to enhancement 
of the cross sections for the W^W~ and ZZ final states, as generically expected for a 
higgsino-like LSP However, it is difficult to see the effect of CP violation, since 
the deviation of the solid line from the dashed one is very small in Fig. ^1 In the case 
of a bino-like LSP as in Fig.|3J the neutralino pair annihilation cross section is smaller 
than that for the case of a mixed LSP. Therefore the supersymmetric contribution to 
the positron flux is completely invisible. 
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^ In the p resent analysis, we have not included coannihilation effects [35 

in the calculation of the relic density. They are crucial in the case that the 
next-lightest superparticle is nearly degenerated with the LSP. Since this situation 
happens when the LSP is higgsino-like, our results on the relic density for a higgsino- 
like LSP may be altered if coannihilations are taken into account. The investigation 
of the effect of CP violation on the coannihilations is left for future work. 
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VII. CONCLUSIONS 



We have studied pair annihilations of the neutrahno dark matter in the minimal 
super symmetric standard model with CP violation. We have considered the case 
that the higgsino mass and the trilinear scalar couplings have CP-violating phases 
of order unity, taking a scenario that the scalar fermions in the first two generations 
are much heavier than those in the third generation to avoid a severe constraint from 
experimental limits on electric dipole moments. It has been found that, when the 
lightest neutralino is bino-like, the cross section of the process XX ~^ W^W~ and 
XX ~^ H2H2 for nonrelativistic neutralinos can be significantly enhanced by 6^, the 
phase of the higgsino mass. It follows that the relic density of the neutralino can be 
considerably suppressed by this effect. However, even this enhancement is not enough 
to make bino-like dark matter consistent with the cosmological constraint ()30p. We 
have also discussed the effect of CP violation on the positron flux from neutralino pair 
annihilations in the galactic halo. It has been shown that the effect of CP violation 
is difficult to see in the positron flux. 
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FIG. 1: Feynman diagrams for xx ~^ ff- the s-channel neutral Higgs boson exchange, the 
s-channel Z-boson exchange, and the t- and u-channel sfermion exchange. Note that the 
u-channel diagram is not shown here. 




FIG. 2: Feynman diagrams for xX ~^ W'^W : the s-channel neutral Higgs boson exchange, 
the s-channel Z^-boson exchange, and the t- and u-channel chargino exchange. 




FIG. 3: Feynman diagrams for xX ~^ -^2-^2' s-channel neutral Higgs boson exchange 
and the t- and u-channel neutralino exchange. 



25 



mo = 1 TeV. = 500 GeV, |A| = 1 TeV 

tan/S = 5. Mg = 300 GeV, = 400 GeV, 6^ = 
0.004 




_Q 001 ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' 

0.5 1 1.5 2 

FIG. 4: The variation of the quantities and defined in eq. (|2U() with the CP 
violating phase 0^. The relevant parameters are taken as mo = = 1 TeV, mA = 500 GeV, 
tan/3 = 5, M2 = 300 GeV, |^| = 400 GeV and Oa = 0. The dashed and solid curves 
correspond to the result for C^i and D^^, respectively. In this case, the LSP is bino-like. 
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mo = 1 TeV. = 500 GeV, |A| = 1 TeV 
tan/S = 5. Mg = 300 GeV, = 400 GeV. 6^ = 
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FIG. 5: The absolute value of the coupling constants C^^^° (r = 1, 2, 3) normalized 
mw as a function of 0^ for the same choice of parameters as Fig. 0] 
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mo = 1 TeV. = 500 GeV, |A| = 1 TeV 




FIG. 6: The imaginary part of the pure imaginary couphng constants (r = 1, 2, 3) 

as a function of 0^ for the same choice of parameters as Fig. H) 
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FIG. 7: The cross section times relative velocity av for v = 10~^ versus for the same 
choice of parameters as Fig. 0] The solid, long dash-dot, short dashed, short dash-dot and 
short dash-long dash lines correspond to the contributions from W^W" , ZZ, bb, t'^t~ 
and H2H2 final states, respectively. The bold solid line represents the sum of all the 
contributions. 
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nio = 1 TeV. = 500 GeV, |A| = 1 TeV 
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FIG. 8: The relic density n^h^ versus 6^ for the same choice of parameters as Fig. 0] In 
the shaded region, the rehc density is consistent with the WMAP 2a bound. 
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FIG. 9: The cross section times relative velocity av for v = W ^ versus 6a for the same 
choice of parameters as Fig. |3]but 6^ = 0. 
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FIG. 10: The cross section times relative velocity av for v = 10^^ versus 6^ for M2 = 
300 GeV and |^| = 200 GeV. The values of the other parameters are the same as Fig. 0] 
The solid, long dash-dot, short dashed, long dashed, short dash-dot, short dash-long dash 
and dotted lines correspond to the contributions from W^W^ , ZZ, bb, tt, t~^t~, H2H2 
and ZH2 final states, respectively. The LSP in this case is a mixture of the bino and the 
higgsinos. 
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rrio = 1 TeV, = 500 GeV. |A| = 1 TeV 

tan/? = 5. M„ = 300 GeV. |m| = 200 GeV, 0=0 

IQO 1^-^ ^^-n ^ ^ ^ ^ ; 
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FIG. 11: The relic density Q^h^ versus 9a for the same choice of parameters as Fig. IIUI 
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riio = 1 TeV. = 500 GeV, |A| = 1 TeV 




ICQ 200 300 400 500 600 700 



\^A (GeV) 

FIG. 12: The final states giving the largest contribution to av for v = 10^^ in the (|/i|, 
M2) plane for mo = \A\ = I TeV, niA = 500 GeV, tan^ = 5, 6*^ = and 6^ = 0. The 
regions where the final states W'^W~ , bb, ti and W^H^ give the largest contribution are 
shown in different gray scales. From lighter to darker, each region corresponds to W~^W~, 
bb, ti and W^H^ in this order. In the darkest strip, the relic density is consistent with 
the WMAP 2a result. The white region is excluded by the LEP limit on the chargino mass 
m - > 104 GeV Q and the hghtest Higgs mass m^ro > 113 GeV Q- 
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riio = 1 TeV. = 500 GeV, |A| = 1 TeV 
tan/S = 5, d^ = 7t/2, = ^ 
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FIG. 13: The same as Fig. ^Jbut for 9^ = Tr/2. The regions where the final states W~^W~, 
bb, tt, H^H^ and W^H^ give the largest contribution are shown in different gray scales. 
From lighter to darker, each region corresponds to W^W~ , bb, ti, H2H2 and W^H^ in 
this order. 
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FIG. 14: Variation of the relic density with 9^ normahzed by that for 6^ = 0, 
n^{e^)/ax{9f, = O), as a function of M2 for niQ = \A\ = I TeV, tua = 500 GeV, tan/3 
= 5 and 6a = 0. Varing 6^ in the range < 0^ < it, the rehc density lies between the 
two sohd lines for |^| = 200 GeV. The region between the two dashed hnes represents the 
corresponding result for = 600 GeV. 
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mo = 1 TeV, = 500 GeV, |A| = 1 TeV 

tarijS = 5, Mg = 300 GeV, = 200 GeV, 6^ = 
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FIG. 15: Positron fraction /{e'^ + e^) versus positron energy for the same choice 
of parameters as in Fig. 1101 The dashed and sohd hnes correspond to the results for 6^ = 
and 6^ = 7r/2, respectively. The points with error bars represent the data from HEAT 
measurement 
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mo = 1 TeV, = 500 GeV, |A| = 1 TeV 

tanjg = 5. Mg = 600 GeV, = 200 GeV, 9^ = 

I I I I I 1 1 1 1 I I I I I 1 1 1 1 I I I I I 1 1 1 
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FIG. 16: The same as Fig.^Jbut for M2 = 600 GeV. The LSP in this case is higgsino-hke. 
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